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ABSTRACT 


In  this  paper  we  are  interested  in  the  problem  of  convergence  of  the  eigenvalues  and  eigenfunctions  of  the 
Laplacian  with  Went  zell-Robin  boundary  condition  to  the  eigenvalues  and  eigenfunctions  of  the  Laplacien  with  Dirichlet 
boundary  condition  when  the  Robin  parameter  tends  to  infinity.  We  show  in  particular  that  the  convergence  of 
eigenfunctions  is  better  than  the  usual  internal  convergence  in  and  boundary  convergence  in  L2(F). 
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1.  INTRODUCTION 


Let  £2  be  a  sufficiently  smooth  domain  (at  least  C2)  and  let  rD  —  Ui<j<ND^Dj 
Ui<j<NRffij  and  rWR  =  Ui<j<NwRf/vflj  respectively  denote  regions  where 

Dirichlet  BC,  Neumann  BC,  Robin  BC  and  Went  zell-Robin  BC  are  imposed. 


^Ne  —  Ui<j<NN^Vej> 


Figure  1:  An  Example  of  an  Open  Domain  with  Different  Conditions  on  Boundary 
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The  question  we  want  to  answer  in  this  article  is: 
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«  For  a  Fixed  Domain 

How  Do  Eigenvalues  and  Eigen  Functions  of  the  Went  zell-Robin  -  Laplacian  Depend  on  the  Parameters  of 
the  Boundary  Conditions  ?  » 

The  Robin  boundary  condition  (also  known  as  Fourier  boundary  condition  or  of  the  third  kind)  and  Went 
zell-Robin  boundary  condition  are  called  artificial  conditions  and  play  important  roles  in  many  physical  phenomena 
governed  by  partial  differential  equations.  We  can  cite  as  examples:  thin  layers  problems,  electromagnetic,  acoustics, 
vibration,  heat  problems. etc.  Went  zell  condition,  which  was  introduced  in  [1],  can  be  understood  as  an  additional  potential 
energy  from  the  boundary  5  fga\Vgnu\2  which  adds  to  the  internal  potential  energy  through  the  outflow.  In  the  dynamic 

case  it  expresses  an  energy  supplement  of  both  potential  and  kinetic  kind:  (Jan|9tu|2  +  5  Jan|Vgnu|2). 
The  parameters  a  and  /?  appearing  in  the  Wentzell-Robin  condition  are  either  operators  or  simply  constants  imposed  by  the 
physical  nature  and/or  by  the  geometric  nature  of  the  problem  considered. 

There  is  a  large  literature  on  the  subject,  we  cannot  make  a  complete  list  in  this  paper,  but  for  some  aspects  of  the 
physical  interpretation  as  well  as  mathematical  and  numerical  study  of  such  conditions  the  following  articles  and  their 
references  can  be  consulted:  [2],  [3], [4], [5], [6],  [7],  [8],  [9]. 

2.  PRELIMINARIES  AND  NOTATIONS 

Let  us  consider  the  Dirichlet  eigen  value  problem  of  the  Laplacian  in  a  C 2 -bounded  domain  ft  c  R2: 
— Au  —  Au  in  H  (5) 
u  =  0  on  dfl  (6) 
There  exist  a  sequence  (Ak,uk)k-1<J3  solutions  of  the  system  and  such  that: 

0  <  A-l  <  A2  <  A3  ...  <Ak  <  ■••  <  ■••  (7) 
Counted  with  their  multiplicities,  and  (  uk)k=laD  is  ortho  normal  in  L2(£2)  satisfying: 

Ak  =  TnaxA^Hi^minvSAv^0codim(A-)=k_1R(k,D,v~)  (8) 
Where  R(A,  D,  v)  denotes  the  classical  Rayleigh  quotient  i.e.  : 

R^,D,v)  =  ^  (9) 
Consider,  on  the  other,  the  following  Went  zell-Robin  eigenvalues  problem: 

-Au  =  AuinH  (10) 
au  +  —  —  SAau  —  0  on  d!2  (11) 
There  exist  a  sequence  (Ak(a,  8),uk(a,  5))fc=l03  solutions  of  the  system  and  such  that: 

0  <21(a,5)  <A2(a,8)  <A3(a,S)...  <Ak(a,S)  <  ■■■  <  ■■■  (12) 
And  we  can  again  make  the  family  (  uk{a,  5))k=lcc  orthonormal  in  L2(0),  and  such  that  the  eigenvalues  satisfy: 
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Ak(cc,8~)  —  maxAaK(a)  minv£A iV*o,codim(A)=k-i  ct,  8,  v)  (13) 
Where 

H(£l)  =  jv  £  H^ifl);  U\r.  £  HliTi)  V  1  <  i  <  Nv;  u\r.  =  0, 1  <j  <  JV0]  (14) 
and  R(A,  a,  §,v)  denotes  the  corresponding  Rayleigh  quotient  i.e.  : 

R(A,  a,  5,  v)  =  v—  v-   (15) 

We  give  in  Section  3  a  an  answer  to  the  question  posed  in  the  introduction  which  can  be  seen  both  qualitative  and 
quantitative.  In  Section  4  we  expose  the  proof  of  the  main  result  in  several  steps.  Finally  we  conclude,  in  section  5,  this 
paper  with  some  remarks,  perspectives  and  some  open  problems. 

3.  MAIN  RESULTS 

We  give  below  the  first  result  concerning  the  dependence,  on  parameters  a  and  8,  of  the  eigenvalues  and 
eigenfunctions  of  the  Laplacian  with  Went  zell-Robin  boundary  condition  and  the  main  result  of  convergence  of  the 
eigenfunctions  when  Robin  parameter  tends  to  infinity: 

Theorem  1 

(a)  Vfe  >  1,  the  Went  zell-Robin  eigenvalue  Xk{a,  5)  is  continuous  converge  to  2.k(D)  when  a  -*  +co  and 

satisfy  : 

Va1  >  a2;  Afc(a1(  8)  >  Ak(a2, 8)  (16) 
V81>82;Xk(cc,8J>Xk(a,82)  (17) 
Va  >  0;Ak(D)  >  Ak(a,8)  (18) 

(b)  V/c  >  1,  the  Wentzell-Robin  eigenfunction  uk(a,  8)  converge  to  uk(D)  in  the  following  sense  : 
uk(a,8)  ->  uk(D)  intf2(Q)  (19) 
uk  (a,  8)  ->  0  in  H1  (TD  U  rfi  U  l»  (20) 
and 

?»^^Minl2(rJ|Ur,)  (2i) 

Recall  that  we  had  already  shown  in  [8]  and  [9]  the  convergence  of  eigen  functions  in  H2(£i)  m  the  case  where 
the  boundary  condition  is  a  pure  Robin  BC  to  the  eigen  functions  of  Dirichlet  when  the  parameter  a  tends  to  infinty  and  we 
had  shown  the  convergence  to  the  eigen  functions  with  Neumann  BC  when  the  parameter  a  tends  to  zero  but  with  the 
restriction  to  the  case  of  convex  domains.  In  this  context  Theorem  1  is  an  extension,  in  the  Dirichlet  case,  of  this  above 
result,  however  here  we  do  not  treat  the  case  of  convergence  to  the  Neumann  problem. 

Note  also  that  in  another  work  being  published  we  got  finer  estimates  for  the  principal  eigen  value  and  its 
asymptotic  expansion  as  a  function  of  a. 
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4.  PROOF  OF  THE  MAIN  RESULTS 

The  proof  is  divided  into  several  steps 
4.1.  The  Eigenvalues  as  a  Function  of  Parameters  Went  Zell-Robin 

There  is  no  difficulty  in  noting  that 

R(A,  a^S,  v)  <  R(A,a2,8,v)  Va±  >  a2  >  0,VS  >  0 

and 

R(A, a,8i,v)  <  R(A,a,82,v)  VS1  >  S2  >  0,Va  >  0 

On  the  other  hand,  for  any  fixed  v,  the  Rayleigh  quotient  is  clearly  a  continuous  function  with  respect  to 
parameters  a  and  5. 

Of  course   both   the   continuity   inequalities   remain   true   through   the   introduction   of  the  minimum: 
min  and  then  the  maximum:    max  . 

v£A,v*0,codim(A)  =  k-l  AcH{0) 

On  the  other  hand  we  have 

Jn|vw|*+gJ  M'+aJ  |vrwi' 
R  (A,  a,  8,  w)  =   > 


R(A,w)  =  ^Vw£//0H/2) 

Jqiwi 


And  since  the  space  //q(/2)  is  included  in  J-Ciii)  the  inequality  Va  >  Q;Xk(D )  >  Xk{a,S)  is  immediate. 

This  completes  the  proof  of  the  first  part  of  theorem  1 . 

4.2.  Convergence  of  Went  Zell-Robin  Eigen  Functions 

The  proof  is  also  done  in  several  steps:  We  first  begin  by  establishing  identities  Rellich  type  for  eigenfunctions 
inspired  by  a  technique  usual  multipliers  control  theory  (Lemma  1  and  Lemma  2),  then  in  the  last  step  we  analyze  the 
passages  to  the  limit. 

Stepl 

Lemma  1 

For  any  real  parameters  a  and  S  and  any  eigen  function  uk{a,  8)  we  have: 
In{\dxxuk(a,S)\2  +  \dyyuk{a,S)\2  +  2\dxyuk{a,S)\2^ 
+  frRurvUdruk(a' 5~)\2  +  \dnuk(a,8)\2}{rxdrTy  -  TydrTx) 
+2a  jrRurv\druk(a,S)\2  +  2  fry\VrdTuk(a,  S)\2  = 

(^)2  In\(uk(a,8))\2  (22) 
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Proof  of  Lemma  1 


We  multiply  the  Went  zel-Robin  eigen  value  equation  by  —  huk{a,  8)  and  integrate.  We  then  obtain: 


L\Auk(cc,8)\2  = 


Ak(a,S)  jn\Vuk(a,S)\2  +  SXk(a,S)  JajVruk(a,5)|2  +  aXk(a,8)  Jn  |(uk(a,5))2 


(23) 


Moreover  one  has 


I4(a.  5)l2/n|"w 


=  L  |Aufc(a,S)|z 


(24) 


and 


L|Aufc(a,S)|2  = 


In\dxxuk(a,  8)\2  +  \dyyuk(a,8)\  +  l(dxxuk(a,8))  (dyyuk(a,8)) 


(25) 


Let  us  develop  the  last  term  as  follows  : 


In2(dxxuk(a,S))  (dyyuk(a,8)j  =  fn\dxyuk(a,S)\2  + 
Saa{.dxuk{a,8)rix)(dyyuk{al8is)  -  {dxuk{a,8)){dxyuk{a,8)r)y) 


(26) 


Define  /  as  the  boundary  integral  in  the  above  identity.  To  do  the  calculations  on  the  boundary  we  use  an  open 
cover  (rj)1=1  N  such  that  Fi  =  {(x,  y);  y  =  gi(x)}  and  where  g;  represents  denote  local  chart  of  class  C2,  and  such  that  we 


Integrations   by  parts   are,   then,   performed   classically   by   straightening   each  part   of  the  boundary. 
We  then  get:    /  = 


k  =  JMTxTydTuk(a,S)dT2Ufc(a,5);/2  =  Jgn  Tx\dT  uk(a,  8)\2dTTy; 
7s  =  SdnTy2  dnuk{a,8)dT2  uk{a,8);l4  =  fgnTydTTydnuk(a,8)  dTuk(a,8); 
7s  =  IgnTx2  dT  uk(a,S)dzdnuk(a,S);  I6  =  fgnTxdTTxdTuk(a,8)dnuk(a,8); 
ll  =  IgnTxdTTxdTuk(a,S)dnuk(a,S)  and/8  =  fgn  ry  dTTx(dnuk(a,  S))2 
Then 

7i  =  ~  hnTxTydTuk{oc,8)dT2Uk{a,8) 
www.iaset.us  editor  @iaset.us 
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=  ~~  j/S/3T*TydT|(dT"te(a,5))|2 

=  ~fgnKdTuk(a,8))\2(TxdT(Ty)  +  TydT(Tx))  (28) 
h  =  Jdn\dTuk(a,8)\2(TxdTTy)  (29) 
k  +  h=  lSdn\drUk(a,8)\2(rydt(Tx)  -  TxdT(ry))  (30) 
h  +k  =  lfgn\dnuk(a,8)\2(TydT(jx)  -TxdT(jy))  (31) 
7s  =  !anTy2dnUk(oc,8)dT2Uk(a,8) 

=  -  jgnTy2dT(dnuk(a,8))dTTy2  (32) 
7s  =  ~  JdnTx2drUk(a,8)dnuk(a,8)  (33) 
One  obtain 
Is  +  h  = 

jgndTuk(a,8)dT((dnuk(a,  8)){-xx2  -  xy2)  -  ^gn{dnuk{a,8)){dTuk{a,8))dT{Ty2) 

=  ~  Jdndruk(a'S)dr<idnuk(a,8))  -  fgn(dnuk(a,8))druk(a,8)dr(jy2)  (34) 
h  +  h  =  lJgndrUk(a,8)(dnuk(a,8))(dTTy2  -  dTxx2)  (35) 
h  +  /S  +  h  +  h  =  B 

=  ~  IgndTuk(a,8)dT(dnuk(a,8))  -  ^jgndTuk(a,8)(dnuk(a,8))dT(Ty2  -t2) 

=  -  jgndTuk(a,8)dT(dnuk(a,8))  (36) 
h  +  /2  + 17  +  I8  = 

\igni.i.dTuk{a,8))2  -  (dnuk{a,8))2)(TydTTx-TxdTTy)  (37) 
And  finally  we  have 
B  =  -  JdndTuk(a:S)dr(dnuk(a,8)) 
=  ~  ia/2  ar"k(a.  8)dr(Aruk(a,  8)  -  auk(a,  8)) 
=  ~  hn  dTukioc,  8)dzAruk(a,  8)  +  a  Jgn  dxuk{a,  8)dTuk(a,  8) 
=  ~  Jdn\ArdrUk(a,8)  |2  +  a/M(aTufc(a,5))2 
=  Jdn\ArdrUk(a,8)  |2  +  a/M(aTuk(a,5))2 

=  JM|aT2Uk(a,5)|2+a/M(aTufc(a,5))2  (38) 


www.iaset.us 


editor  @iaset.us 


On  the  Convergence  of  Eigenvalues  and  Eigen  Functions  of  the  Laplacian  with  Went  Zell-Robin  Boundary  Condition 


17 


Step  2 

Lemma  2:  Let  (<7k)i<fc<n  be  a  vector  field  in  Wl(£i)  such  that  Xi<k<n  <7knk  =  1  on  dQ.,  then  we  have  the  following 
identities 

l^vry^r^'S)]2  -  \dnuk{a,8)\2}  = 
-  L\Vuk{a,  5)|2  — -  L  9uk(a'5)  dqi  dUkia'S) 


dx[      dx[  dxi 


2 


T1  'a  I  C«*  <«'  +  K  C".  ^))2 1  (39) 

and 

-\STR„Tv\dnuk(D)\2  = 

i/jVufcW|^-/^g^-^/n|(ufcW)2|g  (40) 

Sketch  of  the  Proof  of  Lemma  2 

The  identities  of  Lemma  2  are  easily  obtained  by  multiplying  successively  the  kth  eigenvalues  Laplace-Went 
zell-Robin  problem  by  uk  (a,  5)  and  the  kth  eigenvalues  Laplace -Dirichlet  problem  by  uk  (D)  and  in  tegrating  by  parts. 

Step  3 

Note  first  that  there  is  no  particular  difficulty  to  show  for  every  fixed  and  positive  5  that  the  increasing  and 
bounded  sequence  Xk(a,  8)  converges  to  2.k(D).  The  fact  that  the  sequence  uk(a,S)  is  bounded  in  Hx(fi)  then  deduce 
immediately.  So  it  is  easy  to  identify  the  limit  and  to  prove  that  uk{a,  8)  converges  in  H1(££)  to  ufc(£>)  when  a  tends  to 
infinity.  To  be  convinced  we  acn  use,  for  example,  the  following  identity  : 

/Q|v(uk(a,5)  -  uk(D))\2  =  Xk(a,8)fa\uk(a,8)  -  uk(D)\2  + 

(ak(a,8)  -  ak(D))  jQ\uk(D)\2  ~  /9nu^a'5)^  <41> 

Now  observe  that  in  the  case  where  the  quantity  (rxdrTy  —  TydrTx  )  is  positive  then  the  convergence  uk(a,  8)  -* 
0  in  is  immediate.  Note  that  this  quantity  is  positive  especially  when  the  domain  is  convex  [9]. 

In  the  general  case,  we  know  from  the  first  identity  of  Lemma  2  that  there  exists  a  constant  C  independent  of  the 
parameters  of  the  problem  and  such  that: 

SrRurv\dnUk(a,S)\2  <  c{frRurv\dTuk(a,S)\2  +  /JVuk(a, S)\2}  (42) 

Using  both  Lemma  1,  the  above  estimate  and  the  convergence  of  the  sequence  (uk(a,  5))fc  in  H1^)  we  we  are 
able  to  prove  the  convergence  of  the  trace  uk(a,  5)\gn  in  H1(dii)  to  zero. 

To  prove  the  convergence  of  the  normal  derivative,  we  first  begin  by  identifying  its  weak  limit  using  adequate 
raising  "R  from  //1/2(90)  into  HX(Q)  such  that  we  have: 

jgndnuk(a,8)v  =  IQVuk(a,S)V(jl(v))-Ak(a,8)fnuk(a,8)(jl(v))  (43) 
www.iaset.us  editor  @iaset.us 
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w  £//1/2(an) 

This,  of  course,  implies  the  convergence  of  dnuk(a,  8)  in  H1^2  (3£X)  to  a  certain  /  satisfying: 
Jan/17  =  JnVufc(fl)V(52(v))  -  Ak(D)  fnuk(D)(X(v)) 

=  ;an3„uk(Z))i;  (44) 
Vv  £  tf1/2(dQ) 

From  the  estimate  (42)  we  deduce  that  the  convergence  of  the  normal  derivative  is  actually  weak  in  L2(dQ). 
The  convergence  of  the  norm  is  then  obtained  by  simple  passage  to  the  limit  in  the  first  identity  of  Lemma  2.  It  is,  then, 
sufficient  to  identify  term  by  term  the  result  of  the  passage  to  the  limit  obtained  with  the  second  identity  of  Lemma  2  to 
conclude  that  convergence  of  the  normal  derivative  is  strong. 

Finally,  the  convergence  H2  (£2)  becomes  immediately  from  the  identity  of  Lemma  1  and  the  convergence  result 
established  above. 

5.  CONCLUSIONS 

In  this  article,  we  established  the  convergence  result  of  eigenvalues  and  eigenfunctions  of  Laplacien  operator  with 
Went  zell-  Robin  boundary  condition  when  the  Robin  parameter  tends  to  infinity. 

It  may  be  interesting  to  be  able  to  solve  the  problem  of  convergence  to  the  Neumann  problem  when  Robin 
parameter  tends  to  zero.  It  would  also  be  interesting  to  obtain  results  of  asymptotic  behavior  both  in  relation  to  the  setting 
of  the  boundary  conditions  in  relation  to  the  variation  of  geometry. 

Another  interesting  problem  would  be  to  examine  the  convergence  of  Went  zell  control  system  to  the  Dirichlet 
control  system  (respectively  Neumann)  when  the  parameter  a  tends  to  infinity  (respectively  to  zero)  and  to  extend  the 
results  obtained  [9]  in  the  case  of  Robin-Fourier  control. 
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